Non-equilibrium transitions in fully frustrated Josephson junction arrays 
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We study the effect of thermal fluctuations in a fully frustrated Josephson junction array driven by 
a current / larger than the apparent critical current I C (T). We calculate numerically the behavior of 
the chiral order parameter of Z2 symmetry and the transverse helicity modulus (related to the U(l) 
symmetry) as a function of temperature. We find that the Z2 transition occurs at a temperature 
Tz 2 (I) which is lower than the temperature Tt/(r) (-0 for the U(l) transition. Both transitions could 
be observed experimentally from measurements of the longitudinal and transverse voltages. 



PACS numbers: 74.50+r, 74.60. Ge, 74.60.Ec 

The study of non-equilibrium steady states of driven 
many-degrees-of-freedom systems are of importance in 
condensed matter physics . Examples of this prob- 
lem are the dynamics of vortices in type II superconduc- 
tors jl]-[§ and charge density waves In two dimen- 
sions, Josephson junction arrays (JJA) are a well con- 
trolled system fl where this issue can be investigated 
HQ . In the presence of a magnetic field such that there 
is a half flux quantum per plaquette, / = iZa 2 /$o = 1/2, 
the JJA corresponds to the fully frustrated XY (FFXY) 
model ]^-pd| . The ground state is a "checkerboard" vor- 
tex lattice, in which a vortex sits in every other site of 
an square grid || . There are two types of competing or- 
der and broken symmetries: the discrete Z2 symmetry 
of the ground state of the vortex lattice, with an asso- 
ciated chiral (Ising-like) order parameter, and the con- 
tinuous U{1) symmetry associated with superconducting 
phase coherence. The critical behavior of this system has 
been the subject of several experimental jl(| and the- 
oretical [pp|Jll|-p^| studies. There are a Z2 transition 
(Ising-like) and a U(l) transition (Kostcrlitz-Thouless- 
like) with critical temperatures Tz 2 > Tuns. There is 
a controversy about these temperatures being extremely 
close @ or equal QjjJQ . In the light of this, it is 
worth studying the possibility of non-equilibrium Zi and 
t/(l) transitions at large driving currents. Also, the dy- 
namical transitions in driven systems studied up to now 
involve continous (translational or gauge) sym- 
metries, and therefore it is interesting to study a system 
with a discrete symmetry. Previously, we have found dy- 
namical transitions of the vortex lattice in a JJA with a 
field density of / = 1/25 0: for large currents I there 
is a melting transition of the moving vortex lattice at a 
temperature higher than the transverse superconducting 
transition: Tm{I) > Tr/m(J). Interestingly, here we find 
that the opposite case occurs in the driven FFXY: the 
order of the "checkerboard" vortex lattice is destroyed at 
a much lower temperature than the transverse supercon- 
ducting coherence, Tz 2 {I) < Tu(i){I). 

The hamiltonian of the frustrated XY model is: 

n = E ^ cos ^ n + - e ^ A»] , (i) 



where Iq is the critical current of the junction between 
the sites n and n + \i in a square lattice [n = (n x , %), 
fx = x, y] , Rn is the normal state resistance and 0^(n) = 
0(n + /1) - 0(n) - A^n) = A M 0(n) - A^n) is the gauge 
invariant phase difference with A^(n) = |2L j^ + ^ a A-dl. 
In the presence of an external magnetic field H we have 
A^xA^n) = A x (n)-A x {n + y) + A y (n + x.)-A y (n) = 
27r/) / = Ha 2 /$o and a is the array lattice spacing. For a 
fully frustrated JJA we have / = 1/2. Here, we will take 
periodic boundary conditions (p.b.c) in both directions 
in the presence of an external current I ex t in arrays with 
L x L junctions ||. The vector potential is taken as 
-A M (n, t) = A®(n) — a^(t) where in the Landau gauge 
A x (n) = -2-nfn y , A y {n) = and a^(t) will allow for 
total voltage fluctuations. With this gauge the p.b.c. for 
the phases are: 9(n x + L, n y ) = 9(n x , n y ) and 9{n x , n y + 
L) = 0(n x ,n y ) — 2-7T fLn x . The current flowing in the 
junction between two superconducting islands in a JJA 
is modeled as the sum of the Josephson and the normal 
currents |,|,|ll|,0: 



J^(n) = I o sin0 M (n) 



$0 dfl M (n) 
2itcRn dt 



^(M) (2) 



where the thermal noise fluctuations ry M have correlations 
(»7 /1 (n,t)v(n / > 0) = f^ M ' f mn'*(t-t')- The condi- 
tion of a current flowing in the y- direction: zZnJpi 11 ') = 

After 



determines the dynamics of a M (t) 

A M • 7 M (n) = 
we obtain the full RSJ- 



lextL Sfj 

considering local conservation of current, 
E^(n) -Un-fi) = 0, 
Langevin dynamical equations of the driven XY model 
as in (^J^|. We normalize currents by Iq, time by 
tj = 2ttcRnIq/$q and temperature by I^o/2-KkB- We 
solve the dynamical equations with time step At = 
0.001 — O.Itj and integration times 10000tj after a tran- 
sient of 5000tj. 

We study the fully frustrated JJA for system sizes of 
L x L junctions, with L = 8, 16, 24, 32, 48, 64, 128. In 
the absence of external currents, we find an equilibrium 
phase transition at T c = 0.45 which, within a resolution 
of AT = 0.005, corresponds to a simultaneous (or very 
close) breaking of the U(l) and the Z2 symmetries. Here 
we will analyze the possible occurrence of these transi- 
tions as a function of temperature when the JJA is driven 
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by currents well above the zero temperature critical cur- 
rent I > J c0 = (V2 - 1)7 w 0.414J . 



from an ordered state at T = 0, see Fig. la. 
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FIG. 1. Breaking of the (7(1) symmetry for a large current: 
I = 0.9, I > I c {0), system size 64 x 64. a) Helicity Modulus 
T x vs temperature T (• increasing T, O decreasing T). Inset: 
size effect for L = 32, 48, 64, 128. b) Transverse voltage for a 
small transverse current, Itr = 0.1, vs T. c) Vortex-antivortex 
pairs density, n va vs. T. 

U(l) symmetry and transverse superconductivity. In 
the driven JJA superconducting coherence can only be 
defined in the direction transverse to the bias cur- 
rent |7|]l8f|. We calculate the transverse helicity modu- 
lus = ^(E„cos^(n)) - ^{(EnSin^n)] 2 ) - 
(En s ' n ^( n )]) 2 }- P n or der to calculate the helicity 
modulus along x, we enforce strict periodicity in 9 by 
fixing a x (t) = 0]. We find that T x is finite at low T and 
vanishes at a temperature (-0-In Fig. la we show the 
behavior of T X (T) for a current I = 0.9 in a 64 x 64 JJA. 
The inset of Fig. la shows T, for sizes L = 32, 48, 64, 128, 
we see that a transition temperature can be defined in- 
dependently of lattice size. This transition is reversible: 
we obtain the same behavior when decreasing T from a 
random configuration at T = 1 and when increasing T 
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FIG. 2. Breaking of the Z2 symmetry for a large current: 
/ = 0.9, I > I c (0), system size 128 x 128, results for increas- 
ing T (•) and decresing T (O). a) Chiral order parameter \ 
vs T and one-dimensional order parameter \ x vs. T. Inset: 
size effect for X for L = 8, 16, 32, 48, 64, 128. b) Longitudinal 
voltage V vs. T. Insets: snapshots of the staggered mag- 
netization m s (n,t): ordered state for T — 0.035 (warming 
up), high temperature disordered state, T = 0.15, and low 
temperature state with /^-domain walls, T = 0.0025 (cool- 
ing down), c) Domain wall lengths C x and C y vs. T. Inset: 
domain anisotropy C x /C y vs. T. 

Transverse superconductivity can be measured when a 
small current It r is applied perpendicular to the driving 
current: we find a vanishingly small transverse voltage 
Vt r below Tiy(i) (-0 1 as we found before in M for / = 1 /25. 
We obtain the voltage in the /i-direction as the time av- 
erage V M = (da^(t) I dt) (normalized by RnIq)', longitudi- 
nal voltage is V = V y and transverse voltage is Vj r = V x . 
In Fig. lb we see that the transverse resistance Vt r /hr is 
negligibly small for T < Tf/m and starts to rise near the 
transition. The equilibrium U(l) transition (at / = 0, 
1 = 0, Kosterlitz-Thouless) is characterized by the un- 
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binding of vortex- antivortex pairs above T c . We calculate 
the density n va of vortex-antivortex excitations above 
checkerboard vortex configuration as 2n va = (|fe(n)|) — /, 
where the vorticity at the plaquette n (associated to the 
site n) is fc(fi) = — A p x nint[0 M (n)/27r]. We see in Fig.lc 



that 



rises near Ti 



(7(1) 



Moreover, the transverse resis- 



tivity above T^m is Vtr/Itr n M . 

Zi symmetry. Since the ground state is a checkerboard 
pattern of vortices, we define the "staggered magneti- 
zation" as m s (h,t) = (-l)"*+ ,l w [2b(n x , n y , t) — 1] and 
m s (t) = (l/L 2 ) £ fi m s (n,i). At T = 0, I = there are 
two degenerate configurations with m s = ±1. Above 
the T = critical current I c q the checkerboard state 
moves as a rigid structure and m s (t) changes sign pe- 
riodically with time. Therefore we define the chiral order 
parameter as x — ( TO sW)- We start the simulation at 
T = with an ordered checkerboard state driven by a 
current / > I c q and then we increase slowly the tem- 
perature. We obtain that the chirality parameter van- 
ishes at a temperature Tz 21 which is smaller than Tjj^, 
as can be seen in Fig. 2a for I = 0.9. This transition 
is confirmed by the size analysis shown in the inset of 
Fig. 2a: for T < Tz 2 the chirality x is independent of 
size, while for T > Tz 2 we see that x ~ l/L 2 . As 
it is shown in Fig. 2b, the longitudinal voltage V has a 
sharp increase at Tz 2 , which could be easily detected ex- 
perimentally. The excitations that characterize the Z2 
transition are domain walls that separate domains with 
different signs of m s . The length of domain walls in the 
direction fi is given by £ M = (2/L 2 ) ^ R (6(n)6(n + i>)), 
with v _L \i. We find that for I > and T > the do- 
mains are anisotropic, with the domain walls being longer 
in the direction perpendicular to the current (C x > C y ) 
and the domain anisotropy C x /C y increasing with /. In 
Fig. 2c we show the dependence of with temperature 
for / = 0.9. At T — there are no domain walls, since 
the initial condition is the checkerboard state, and the 
domain wall length grows with T, showing a sharp in- 
crease at Tz 2 . The domain anisotropy C x /C v is shown 
in the inset of Fig. 2c: it has a clear jump at the transi- 
tion in Tz 2 while for T 3> Tz 2 the domains tend to be 
less anisotropic. When decreasing temperature from a 
random configuration at T = 1, an important number of 
domain walls along the x direction remain frozen below 
Tz 2 - C x tends to a finite value when T — > and the 
domain anisotropy tends to diverge when cooling down. 
This leads to a strong hysteresis in the voltage V at Tz 2 
(see Fig. 2b) since the extra domain walls increase dissi- 
This low T state with frozen-in domain 
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walls is ordered along the x-direction (i.e. perpendicular 
to I) but is disordered along the y direction which gives 
X ~ 0. We define the Z-i order parameter in the x direc- 
tion as Xx = ((!/£) E„ [(V£)E, lx ™sK,%,£)] 2 ) and 
X y is defined analogously. We see in Fig. 2a that, when 
cooling down from high T, Xx vanishes as Xx ~ l/L for 
T > Tz 2 (it has stronger size effects than x) an d becomes 



finite for T < Tz 2 , while Xy ~ for any T. Therefore, 
depending on the history, there are two kinds of high 
current steady states with broken Z2 symmetry at low 
T, examples of which are shown in the inset of Fig. 2b. 
One state has mostly the checkerboard structure (x 7^ 0) 
with few very anisotropic domains. It can be obtained 
experimentally by cooling down at zero drive and then 
increasing /. The other steady state is ordered in the 
direction perpendicular to / (xx 7^ 0, Xy — 0) with sev- 
eral domain walls along the x direction. These domain 
walls move in the direction parallel to / (via the motion 
of vortices perpendicular to /) giving an additional dis- 
sipation. This state can be obtained experimentally by 
cooling down with a fixed /. 
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FIG. 3. Current-voltage hysteresis for T = 0.02 shown for 
(a) voltage vs. / and (b) chiral order parameter vs. /. In- 
creasing current from the checkerboard state (•) and decreas- 
ing current from a random state at large I > Iz 2 (O). 

The two steady states have also different critical cur- 
rents as can be observed in the low T current-voltage 
(IV) characteristics. In Fig. 3a we show the IV curve for 
T = 0.02 and in Fig. 3b the corresponding x vs - I curve. 
When increasing / from the 1 = equilibrium state, we 
find a critical current I C 2(T), which in the limit of T = 
tends to I c o = \/2 — 1 as found analitically and in simula- 
tions with p.b.c |§, pfj| , ^l| . Near I C 2 the order parameter 
X has a minimum and rapidly increases with /. The 
driven state is an ordered state similar to the one shown 
in Fig. 2b. At a higher current Iz 2 there is a sharp drop 
of x which corresponds to the crossing of the Tz 2 (I) line 
(see Fig. 4) and the Z2 order is lost. If we now decrease 
the current either from the disordered state at / > Iz 2 
or from a random initial configuration or from a config- 
uration cooled down at a fixed / > I C 2, we obtain the 
steady state with domain walls along the x direction and 
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X ~ 0, Xx 7^ 0. This state has a higher voltage and pins 
at a lower critical current I C \{T), which has the T = 
limit I c i(T — > 0) = 0.35. It has been shown recently H] 
that open boundary conditions can nucleate domain walls 
leading to the critical current / c i(0) = 0.35 usually found 
in open boundary T = simulations p] , [i7||l9| ] . Also a 
moving state with parallel domain walls (as in the inset 
of Fig. 3b) has been found by Gr0nbech-Jensen et al. ]22j | 
in / = 1/2 J J A with open boundaries and Marino and 
Halsey [p3| have shown that the high current states of 
frustrated JJA can have moving domain walls. We have 
studied the effect of open boundaries in the direction of 
/, in the direction perpendicular to / and in both direc- 
tions. They differ mainly in the T = critical current 
and IV curve, for finite T there are small differences in 
the detailed shape of the hysteresis in critical current. 
In all the cases the two high current steady states are 
observed at finite T with the same history dependence. 
Also, we find that the density of frozen C x domain walls 
depends on cooling rate and decreases with system size. 
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FIG. 4. Current-temperature phase diagram. Tjjci)(I) line 
obtained from T X (T) and V tr (T) (A). T Z . 2 (J) line obtained 
from X (T), V(T) and C x /C y (T) (.). J cl (T) ( O) and J c2 (T) 
(•) are obtained from hysteresis in IV curves as well as from 
hysteresis in T X (T) and x(T) curves. The dashed line corre- 
sponds to the IV curve of Fig.3 (T = 0.02). 

In summary, we have obtained the 
current-temperature phase diagram of the driven fully 
frustated XY model, which is shown in Fig. 4. At high 
currents the breaking of the U(l) and the Zi symmetries 
occurs at well separated temperatures, with Tz 2 < Tj/m- 
The low temperature regime T < Tz 2 (I) has bistability 
with two possible steady states and history dependent 
IV curves. The different transitions could be observed 
experimentally with measurements of the transverse and 
longitudinal voltage. 
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discussions and Fundacion Antorchas and Conicet (Ar- 
gentina) for financial support. 



[1] A. E. Koshelev and V. M. Vinokur, Phys. Rev. Lett. 73, 

3580 (1994); T. Giamarchi and P. Le Doussal, ibid. 76, 

3408 (1996); P. Le Doussal and T. Giamarchi, Phys. Rev. 

B 57, 11356 (1998); S. Scheidl and V. M. Vinokur, ibid. 

57, 13800 (1998); L. Balents, M. C. Marchetti and L. 

Radzihovsky, ibid. 57, 7705 (1998). 
[2] S. Bhattacharya and M. J. Higgins, Phys. Rev. Lett. 70, 

2617 (1993); F. Pardo et al, Nature (London) 396, 348 

(1998). 

[3] K. Moon, R. T. Scalettar and G. Zimanyi, Phys. Rev. 

Lett. 77, 2778 (1996); S. Ryu et al, ibid. 77, 5114 (1996); 

N. Gr0nbech- Jensen, A. R. Bishop and D. Dommguez, 

ibid. 76, 2985 (1996); C. Reichhardt et al, ibid. 78, 2648 

(1997); C. J. Olson, C. Reichhardt and F. Nori, ibid. 

81, 3757 (1998); A. B. Kolton, D. Dommguez and N. 

Gr0nbech- Jensen, ibid. 83, 3061 (1999). 
[4] L. Balents and M. P. A. Fisher, Phys. Rev. Lett. 75, 4270 

(1995). 

[5] For recent reviews see R. S. Newrock, C. J. Lobb, U. 
Geigenmuller and M. Octavio, Solid State Physics 54, 
263 (2000); Lectures on superconducting networks and 
mesoscopic systems, Ed. C. Giovanella and C. J. Lam- 
bert (American Institute of Physics, AIP Conf. Proc. 427, 
New York, 1998). 

[6] D. Dommguez, Phys. Rev. Lett. 72, 3096 (1994) and 82, 
181 (1999). 

[7] V. I. Marconi and D. Dommguez, Phys. Rev. Lett 82, 

4922 (1999), and Phys. Rev. B (2001). 
[8] S. Teitel and C. Jayaprakash, Phys. Rev. B 27, 598 

(1983); Phys. Rev. Lett. 51, 199 (1983). 
[9] B. Berge et al, Phys. Rev. B 34, 3177 (1986). 
[10] B. J. van Wees, H. S. J. van der Zant, and J. E. Mooij, 
Phys. Rev. B 35, 7291 (1987); Ph. Lerch et al., ibid. 41, 
11479 (1990); H. S. J. van der Zant, H. A. Rijken, and J. 
E. Mooij, J. Low. Temp. Phys. 79, 289 (1991). 
[11] K. K. Mon and S. Teitel, Phys. Rev. Lett. 62, 673 (1989). 
[12] J. Lee, J. M. Kosterlitz, and E. Granato, Phys. Rev. B 
43, 11531 (1991); E. Granato et al., Phys. Rev. Lett. 66, 
1090 (1991). 

[13] G. Ramirez-Santiago and J. V. Jose, Phys. Rev. Lett. 68, 

1224 (1992); Phys. Rev. B 49, 9567 (1994). 
[14] P. Olsson, Phys. Rev. Lett. 75, 2758 (1995). 
[15] M. Benakli and E. Granato, Phys. Rev. B 55, 8361 

(1997) . 

[16] E. H. Boubcheur and H. T. Diep, Phys. Rev. B 58, 5163 

(1998) . 

[17] J. S. Chung, K. H. Lee, and D. Stroud, Phys. Rev. B 40, 
6570 (1989); F. Falo et al, ibid. 41, 10983 (1990). 

[18] S. Kim and M. Y. Choi, Phys. Rev. B 48, 322 (1993). 

[19] M. V. Simkin, Phys. Rev. B 57, 7899 (1998). 

[20] M. S. Rzchowski et al, Phys. Rev. B 42, 2041 (1990). 

[21] B. J. Kim and P. Minnhagen, Phys. Rev. B 61, 7017 
(2000). 

[22] N. Gr0nbech- Jensen et al, Phys. Rev. B 46, 11149 
(1992). 

[23] I. F. Marino, T. C. Halsey, Phys. Rev. B 50, 6289 (1994). 



4 



